1/1 


1  AD-H151  513  ON  THE  EFFECTS  OF  DELAMINATION  DAMAGE  IN  FIBRE 

COMPOSITE  LAMINATES(U)  AERONAUTICAL  RESEARCH  LABS 
MELBOURNE  (AUSTRALIA)  R  JONES  ET  AL.  JUN  84 
UNCLASSIFIED  ARL-STRUC-R-403  F/G  11/4  NL 


:OPY  RESOLUTION  TEST  CHART 

NAl  BUREAU  OF  STANDARDS- 1963-A 


■  -m - W •"  ?  a  «*  «* _  «  ,  «fc_  K  .  •*.  a  -  »*,  .  *  *  *  .  *  *  *  m  *  4  *  »  * 

,'*  ,%  ,*»  .*■  %*»  a  -  a  *4  *  *  *  •  *  4  *  >  ’/  ,*  »• 


‘V. 


SIB  FILE  COW  AD-A151  513 


ARL-STMJC-R-403/ 

ARL-MAT-R-116 


AR-003-927 


DEPARTMENT  OF  DEFENCE 

DEFENCE  SCIENCE  AND  TECHNOLOGY  ORGANISATION 
AERONAUTICAL  RESEARCH  LABORATORIES 

MELBOURNE.  VICTORIA 


STRUCTURES  REPORT  403  /  MATERIALS  REPORT  116 


ON  THE  EFFECTS  OF  DELAMINATION 
DAMAGE  IN  FIBRE  COMPOSITE  LAMINATES 


by 

R.  JONES,  J.  PAUL  and  W.  BROUGHTON 


CORY  No 


© 


COMMONWEALTH  OF  AUSTRAIJ 


Ccfrwnoiiv.ea  ih  of  Austro, ,o 

JUNE  1984 


DEPARTMENT  OF  DEFENCE 
AERONAUTICAL  RESEARCH  LABORATORIES 
DEFENCE  SCIENCE  AND  TECHNOLOGY  ORGANISATION 


AR-003-927 


STRUCTURES  REPORT  403  /  MATERIALS  REPORT  116 


ON  THE  EFFECTS  OF  DELAMINATION 
DAMAGE  IN  FIBRE  COMPOSITE  LAMINATES 


by 

R.  JONES,  J.  PAUL  and  W.  BROUGHTON 


SUMMARY 

This  paper  describes  the  results  of  a  numerical  investigation  into  the  effects  of 
delamination  and  impact  damage  on  the  compressive  strength  of  graphite  epoxy  laminates 
and  bonded  metal-to-composite  joints.  For  (he  laminates  considered  it  is  shown  that  as  the 
size  of  the  damage  increases  a  stage  is  reached  after  which  any  further  significant  increase 
in  the  damage  results  in  only  a  relatively  small  decrease  in  the  residual  compressive  strength. 
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1.  INTRODUCTION 


Delamination  damage  in  fibre  composite  materials  may  occur  due  to  a  variety  of  reasons, 
such  as  low  energy  impact  or  manufacturing  defects.  The  presence  of  delamination  damage 
is  of  major  concern  in  the  vicinity  of  bonded  joints  and  in  compressively  loaded  components 
where  damage  may  grow  under  fatigue  loading  by  out-of-plane  distortion. 

An  early  study  into  delamination  growth  arose  from  the  B-l  composite  development 
program  [1].  This  showed  that  delaminations  can  significantly  reduce  the  fatigue  life,  and  the 
residual  compressive  strength,  for  a  compression  dominated  fatigue  load  spectrum.  These 
effects  have  been  confirmed  in  a  series  of  recent  articles  [2,  3,  4,  3]. 

The  present  paper  forms  part  of  a  joint  investigation  into  delamination  damage  currently 
underway  at  both  the  Aeronautical  Research  Laboratories,  Australia  and  the  Royal  Aircraft 
Establishment,  England.  In  this  work  a  three-dimensional  finite  element  analysis  is  performed 
in  order  to  understand  the  mechanisms  involved  in  delamination  growth,  and  subsequent 
failure,  under  compressive  loading.  This  shows  that  as  the  size  of  damage  increases  a  stage  is 
reached  after  which  any  further  increase  in  the  damage  produces  only  a  relatively  small  decrease 
in  the  residual  compressive  strength. 


2.  LAMINATE  CONFIGURATION 

After  examining  the  fracture  surfaces  of  a  quasi-isotropic  laminate,  which  had  failed  under 
a  compression  fatigue  spectrum,  it  was  decided  to  test  a  series  of  laminates  which  contained 
a  teflon  disk  between  the  second  and  third  plies.  The  laminates  were  graphite  epoxy  with  the 
following  lay  up: 

(i)  Type  A:  (0/145/90)38 

(ii)  Type  B:  (0/±45/90)4s 

(iii)  Type  C:  (±45/0/90/0)28 

The  coupons  constructed  from  laminate  type  A  were  76  -2  mm  long,  and  101 -6  mm  wide 
and  contained  a  simulated  delamination  25-4  minx  19-05  mm  between  the  2nd  and  3rd  plies. 
The  coupons  constructed  from  laminate  type  B  were  304-8  mm  long  and  101 -6  mm  wide 
with  a  delamination  38  -1  mm  x  25 -4  mm  in  the  same  location.  These  two  laminates  were  tested 
at  the  RAE  and  details  of  the  experimental  results  can  be  found  in  [20], 

The  coupons  constructed  from  laminate  type  C  were  105  mm  long  and  45  mm  wide  and 
contained  a  variety  of  simulated  delaminations  and  low-energy  impact  damage.  These  coupons 
were  tested  at  ARL  to  determine  the  effect  of  the  damage  on  the  compressive  strength  of  the 
laminate. 

A  detailed  three-dimensional  finite  element  analysis  was  performed  on  each  set  of  coupon 
tests  and  the  results  of  this  investigation  can  be  found  in  the  following  Sections. 


3.  FINITE  ELEMENT  ANALYSIS 

The  finite  element  analysis  of  complex  three-dimensional  delamination  damage  in 
laminated  composites  is  particularly  difficult.  In  general  each  ply  must  be  modelled  separately 
in  order  to  obtain  the  correct  values  for  the  peel  and  the  interlaminar  shear  stresses  around  the 


delamination.  In  addition  a  fine  mesh  is  required  around  the  front  of  the  delamination  in  order 
to  correctly  model  the  stress  singularity.  This  results  in  a  very  large  numerical  model.  In  order 
to  reduce  the  total  number  of  degrees-of-freedom  it  is  tempting  to  use  eight-noded  isoparametric 
bricks.  Indeed  this  approach  was  used  in  [6, 7].  Unfortunately  these  elements  cannot  model 
the  significant  bending  stresses  which  arise  in  the  delamination  problem  [8].  Furthermore, 
since  modelling  is  taking  place  at  the  ply  level  the  elements  have  very  large  aspect  ratios.  This 
gives  rise  to  problems  of  numerical  ill-conditioning.  Again,  using  the  eight-noded  bricks  it  is 
not  possible  to  improve  the  conditioning  of  the  problem.  In  general  these  elements  should  not 
be  used  with  aspect  ratios  greater  than  five  to  one. 

Consequently  under  no  circumstances  should  eight-noded  brick  elements  be  used  to  model 
delamination  problems  in  fibre  composites. 

The  present  investigation  was  done  in  double  precision  using  twenty-noded  isoparametric 
elements  and  a  directionally  reduced  integration  scheme  with  2x2x3  Gaussian  quadrature 
points,  with  3  points  being  used  through  the  thickness  of  the  ply.  This  integration  scheme  is 
described  in  more  detail  in  [9],  However,  as  mentioned  above,  if  each  ply  is  modelled  separately 
the  numerical  model  becomes  excessively  large  and  so  a  new  super-element  was  developed. 

In  this  paper  the  plies  above  the  delamination  are  modelled  separately  as  are  the  two  plies 
below  the  delamination.  The  remaining  plies  are  treated  as  a  super  element  with  the  displace¬ 
ments  varying  quadratically  in  the  local  isoparametric  coordinate  system,  as  in  Reissner  thick- 
plate  theory.  With  this  approach  the  stiffness  matrix  for  each  super-element  can  still  be  written 
in  the  conventional  form,  viz. : 


BDBFdv 


-  i  JJJ 

i=  i  y, 


BDiBPdv 


(D 


Here  D\  is  the  elasticity  matrix  for  the  ith  ply  which  has  volume  V\. 

There  may  be  an  arbitrary  number  of  plies  in  a  single  super-element  and  thus  unlike  classical 
elements  D  varies  throughout  the  thickness  of  the  element.  Details  of  the  Gaussian  quadrature 
required  to  integrate  equation  (1)  correctly  for  an  arbitrary  number  of  plies  is  given  in  [9]  and 
[10].  Such  an  approach  significantly  reduces  the  number  of  nodes  and  elements  required  and 
yet  still  allows  for  an  accurate  calculation  of  the  stress  field.  However,  even  with  this  approach 
the  problem  consisted  of  six  hundred  of  the  twenty-noded  isoparametric  elements.  As  is  now 
standard  practice  the  mid-side  nodes  of  the  elements  surrounding  the  delamination  were  moved 
to  the  quarter  points  in  order  to  simulate  the  required  singularity.  It  must  be  noted  that  when 
using  eight-noded  bricks  as  in  [6,  7]  it  is  not  possible  to  represent  the  required  stress  singularity. 

Having  performed  the  stress  analysis  we  must  now  decide  on  a  criterion  for  assessing  the 
severity  of  the  delamination  damage.  Ideally  we  would  like  to  evaluate  the  energy  release  rate  G. 
This  can  be  readily  done  for  two-dimensional  problems  using  either  the  method  of  virtual  crack 
extension,  as  described  in  [11]  for  the  edge  delamination  problem,  or  in  the  case  of  self-similar 
growth,  by  evaluating  the  J  integral.  For  non  self-similar  growth  /  #  G. 

For  problems  which  are  not  two  dimensional  there  is  no  simple  method  for  evaluating  the 
energy  release  rate  G  without  a  prior  knowledge  of  the  way  in  which  the  delamination  will  grow. 
Furthermore  whilst  a  line  integral  Ji*  has  been  developed  [12],  for  three-dimensional  fracture 
problems,  this  integral  does  not  equal  the  true  local  energy  release  rate.  This  integral  is  defined 
by: 


Ji*  =  limit  ® 
e— o 


( ,  dui  d«i\  , 
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where  the  integration  is  along  a  contour  of  radius  p  normal  to  the  front  of  the  delamination. 
This  integral  coincides  with  the  classical  definition  of  J  for  2-D  problems,  but  for  three- 
dimensional  problems  returns  the  local  energy  required  for  self-similar  growth,  i.e.,  circles  into 
concentric  circles,  and  not  the  true  energy  release  rate.  For  metallic  components  under  mode  1 
fracture  Ji*  is  directly  proportional  to  the  local  stress  intensity  factor  K.  Hence  for  metals  Ji* 
is  a  useful  quantity.  However,  for  our  present  problem  the  growth  is  mixed  mode  and  non  self¬ 
similar  with  the  result  that  Ji*  is  of  questionable  value. 

Alternative  approaches  for  assessing  the  severity  of  cracks  and  holes  in  fibre  composite 
laminates  have  recently  been  developed  [13, 14].  These  methods  are  termed  the  point  and 
average  stress  failure  criteria.  The  point  stress  failure  criterion  assumes  that  failure  will  occur 
when  the  normal  stress,  an,  to  the  crack  (or  hole)  at  a  distance  d0  in  front  of  the  crack  reaches 
the  unnotched  strength  <ju  of  the  laminate,  viz. : 


°n  I,  =  d0  —  CTu  (3) 

The  quantity  d0  is  usually  called  the  damage  zone  size  and  for  graphite  epoxy  laminates  is  typically 
0-9  mm  (=0  038"). 

The  present  paper  uses  the  strain  energy  density  approach  [IS,  16]  to  assess  the  delamination 
damage.  This  approach  may  be  considered  as  an  extension  of  the  point  stress  failure  criterion 
to  allow  for  mixed  mode  failure. 

Define  the  strain  energy  density  in  the  usual  fashion: 


dw/dv  =  ]  cij  (4) 

Then  for  a  two-dimensional  problem  the  strain  energy  density  approach  has  two  basic  hypotheses 
which  apply  for  crack  extension : 


(1)  The  crack  will  grow  in  the  direction  9  —  0o  of  maximum  potential  energy  density 
(viz. :  minimum  strain  energy  density). 


(2)  Failure  occurs  when  the  stress  field  at  a  distance  rc  ahead  of  the  crack  (or  hole) 
in  the  direction  8  =  0O  of  the  minimum  strain  energy  density,  is  such  that: 


{dw/dv)  r,rc  =  (dw/dv) c 


(5) 


Here  ( dw/dv)c  is  the  value,  at  failure  of  the  strain  energy  density  of  the  undamaged  laminate. 
In  this  formulation  rc  plays  a  similar  role  to  d0,  the  damage  zone  size  used  in  the  point  stress 
failure  formulation.  Indeed  as  a  first  approximation  we  can  take  rc  ~  d0.  In  this  approach 
the  strain  energy  density  function  S  defined  by 

§ 

dw/dv  =  -  +  higher  order  terms  (6) 

plays  a  central  role.  For  mode  1  failure  of  orthotropic  material  5  is  proportional  to  the  stress 
intensity  factor  K\.  Indeed  the  above  hypothesis  can  be  readily  expressed  in  terms  of  S;  viz.- 

Failure  occurs  in  the  direction  90  for  which 

8S/89  =  0  and  >  0  (7) 

and  when  the  load  is  such  that 

S/9.9,  =  Sc  (=rc  (dw/dv)e)  (8) 
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For  three-dimensional  damage  we  must  first  locate  at  each  point  along  the  damage  front 
the  local  minimum  of  the  strain  energy  density  function  Smin(=  r  dw/dvm jn).  Failure  then  initiates 
at  the  point  along  the  front  which  has  the  maximum  value  of  Smu». 

For  incremental  growth  each  point  along  the  damage  front  advances  a  distance  r  which  is 
determined  from  the  relationship 

Si  _  Ss  _  Si  /  dw\ 

n  n  '  ‘  ‘  ri  \dv  Jc 

A  more  detailed  description  of  this  hypothesis  including  its  applicability  to  mixed  mode 
crack  growth  can  be  found  in  {15, 161. 


4.  NUMERICAL  RESULTS 


For  the  tests  on  laminates  types  A  and  B  the  finite  element  model  yielded  the  maximum 
values  of  Sam  along  the  lines  AA'  and  BB',  see  Figure  1.  Indeed  the  values*  along  these  lines 
were  relatively  constant  with  the  maximum  value  of  5mio  occurring  at  point  D,  approximately 
3  mm  below  the  centre  line  for  specimens  A  and  B.  Specimen  C,  constructed  from  laminate 
type  C,  had  the  maximum  value  of  Sam  occurring  exactly  at  the  midsides  of  AA'  and  BB'. 

The  values  of  S*  =  (Smin)mmx  for  laminates  type  A  and  B  are  given  in  Table  1  for  a  uniform 
compressive  strain  of  0  004  applied  to  the  ends  of  the  specimen. 


TABLE  1 
Values  of  S* 


Laminate  type 

S*(MPa.mm) 

A 

0  0033 

0  0061 

B 

0  0022 

0  0075 

In  order  to  estimate  the  failure  load  we  need  the  critical  value  of  S,  i.e.  Sc.  One  early  value 
of  Sc  for  an  epoxy  was  given  as  0  055  lb/in  (=  0  0096  MPa.mm),  but  more  recent  work  [16] 
has  found  that  a  value  of  5C  =  0  044  lb/in  (=  0  0077  MPa.mm)  is  more  representative.  This 
value  is  for  a  Modulite  II  5206  graphite  expoy. 

Using  this  value  for  Sc  and  assuming  that  the  coupon  behaviour  remains  linear  elastic 
we  can  now  estimate  the  strain  required  to  cause  failure  which  we  will  denote  as  ff.  This  value 
is  given  in  Table  1. 

It  is  particularly  interesting  to  compare  these  predicted  compressive  failure  strains  with 
those  given  in  [2]  for  a  similar  quasi-isotropic  laminate  with  32  plies,  cf.  laminate  type  B,  and  one 
with  24  plies;  see  Tables  2  and  3  respectively. 


TABLE  2 

Compressive  failure  strain  for  a  32-ply  quasi-isotropic  laminate  [from  (2)] 


Average  Damage  Area  (mm) 

«r 

426 

523 

587 

671 

929 

0  0077 

0  0075 

0  0073 

0  0067 

0  0075 

*  The  values  of  S  were  obtained  via  an  interpolation  procedure  which  is  described  in  the  Appendix. 


TABLE  3 


Compressive  failure  strain  ft  for  a  24-ply  (0/45/0*/— 45/0*/45/0*/— 45/0)s  laminate 

[from  (2)] 


Average  Damage  Area  (mm) 

ft 

458 

574 

645 

703 

955 

0  0060 

0-0064 

0-0057 

0-0057 

0-0060 

As  can  be  seen  from  Table  2  the  failure  strains  predicted  for  laminate  type  B  compare  favourably 
with  those  measured  in  [2]. 

For  the  coupons  constructed  from  laminates  type  C  three  different  sizes  of  delaminations 
were  modelled.  The  delaminated  areas  chosen  were 


(a)  1”  x  1"  (25 -4  mm  x  25 -4  mm) 

( b )  l"xlj"  (25 -4  mm  x  38  •  1  mm) 

(c)  l"x2"  (25 -4 mm x 50 -8 mm) 

In  each  case  the  value  of  ((dw/dv) mln)m»x  occurred  exactly  at  the  centre  of  lines  AA'  and  BB'. 
The  corresponding  values  of  S*  =  (5mm)m»x  for  a  compressive  strain  of  0  004  are  shown  in 
Table  4  along  with  the  predicted  values  for  *r  assuming  as  before  that  Sc  =  0  044  lb/in 
(=  0  0077  MPa.mm). 


TABLE  4 


Delamination  sizes  (mm2) 

645 

967 

1290 

S*  (MPa.mm) 

0  0036 

0-0034 

0-0033 

ft 

0  0057 

0-0060 

0  0061 

Allowing  for  numerical  error  the  value  of  S*,  and  hence  the  failure  strain  «r,  remains  fairly  con¬ 
stant  as  the  delamination  size  increases.  This  phenomenon  can  be  seen  to  occur  in  the  experi¬ 
mental  results  given  in  [2]  and  summarized  in  Tables  2  and  3.  Indeed  this  can  also  be  seen  in 
the  experimental  results  given  in  [4]. 

In  the  coupon  tests  on  laminate  C  a  number  of  6  - 35  mm,  15  8  mm  and  25-4  mm  diameter 
inclusions  were  tested.  These  specimens  failed  at  far-field  strains  of  approximately  5200  ft, 
5100/t  and  4980  /t  respectively.  Whilst  there  is  a  slight  reduction  in  the  far-held  failure  strain 
in  going  from  the  6-35  mm  diameter  inclusion  to  the  25-4  mm  diameter  inclusion,  this  increase 
in  area  resulted  in  only  a  5%  decrease  in  strength. 


In  addition  to  this  set  of  tests  a  second  series  of  tests  was  performed  on  impact-damaged 
specimens.  The  failure  strain  «f  for  these  specimens  is  shown  in  Table  S. 


TABLE  5 


Compressive  strength  of  impact  damaged  specimens 


Damaged  area  (mm2) 

00 

38 

195 

314 

0  0079 

0  0067 

0  0058 

0  0049 

It  thus  appears  that  as  the  size  of  the  damaged  area  increases  the  value  of  S*  (=  (■SmmWx) 
asymptotes  to  a  constant  value.  As  a  result  since  failure  occurs  when  S*  =  Sc  the  residual 
strength  also  asymptotes  to  a  constant  level  as  the  damage  area  increases. 

It  is  important  to  note  that  this  phenomenon  has  also  been  observed  in  the  compressive 
strength  of  composite  laminates  containing  edge  delaminations  [11,  17,  18,  19].  Although  the 
present  paper  has  used  only  small  deformation  theory  in  attempt  to  include  large  deformation 
effects  is  currently  underway.  Indeed  recent  analytical  work*  by  Professor  D.  C.  Stouffer  in 
the  Department  of  Aerospace  Engineering  at  the  University  of  Cincinnatti  has  shown  that, 
even  when  large  deformations  occur,  the  energy  release  rate,  and  hence  ej,  asymptotes  to  a 
constant  value  as  the  size  of  the  delamination  increases. 


5.  DAMAGED  STEP-LAP  JOINTS 

Let  us  now  examine  the  effect  of  delamination  damage  at  a  bonded  step-lap  joint.  In  order 
to  understand  the  mechanisms  involved  a  relatively  simple  joint  configuration  was  considered; 
see  Figure  2.  The  composite  is  a  (±45/02)2s  graphite  epoxy  with  the  delamination  occurring 
at  the  interface  between  the  titanium  step  and  the  0  degree. 

The  moduli  of  the  graphite  epoxy  were  taken  to  be 

G13  =  G12  =  G23  =  5  ’0  GPa, 

£22  =  £33  =  9-5  GPa, 

£11  =  141  GPa, 

V12  =  V13  =  O’ 31, 

V23  =  0-021. 


The  titanium  was  assumed  to  have  £=110  GPa  and  1*  =  0-3. 

As  before  each  ply  was  modelled  separately  using  twenty-noded  isoparametric  bricks  and 
the  region  around  the  crack  front  was  modelled  using  the  fifteen-noded  isoparametric  wedge 
elements  with  the  mid  points  moved  to  the  quarter  points.  A  cross-section  of  the  finite  element 
mesh  can  be  seen  in  Figure  3. 


*  Private  communication. 
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LAYUP  AND  DIMENSIONS 


✓  V 


The  delamination  was  found  to  close  under  applied  compressive  loading.  As  a  result  con¬ 
straints  were  applied  to  the  faces  of  the  delamination  so  as  to  prevent  the  faces  from  crossing 
(i.e.  overlapping).  In  this  case  the  stresses  and  deformations  were  essentially  two-dimensional 
so  that  the  energy  release  rate  G  could  be  calculated  from  the  relationship. 


r  -  — 

~  IB  da 


where  P  is  the  applied  load  B  is  the  width  of  the  specimen  and  8  is  the  movement  of  the  load 
point. 

The  values  of  G  thus  calculated  are  given  in  Table  6  for  various  size  delaminations.  These 
values  correspond  to  an  applied  compressive  strain  of  4000  ji. 


TABLE  6 


Delaminated  area  (mm2) 

31-75 

38-1 

41-27 

44-45 

47-6 

50-8 

G(MPa.mm) 

0-0254 

0-0284 

0-029 

0-0302 

0-0314 

0-0326 

Here  we  see  that  going  from  a  delamination  of  31  -75  mm2  to  50-8  mm2,  a  60%  increase 
in  the  size  of  the  damage,  results  in  an  increase  in  G  of  only  28%.  This  corresponds  to  a  decrease 
in  the  residual  strength  of  only  13%. 


6.  CONCLUSION 

In  this  work  we  have  found  that  for  delamination  damage,  to  both  composite  laminates 
and  metal-to-composite  joints,  as  the  size  of  the  damage  increases  a  stage  is  reached  after  which 
any  further  significant  increase  in  the  damage  produces  only  a  small  decrease  in  the  residual 
compressive  strength. 

The  next  phase  of  this  investigation  involves  detailed  specimen  testing  in  order  to  confirm 
the  numerical  results. 
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